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The article contains general formulas for the production of J = 3/2 resonances by neutrinos 
and antineutrinos. It specializes to the P33(1232) resonance whose form factors are determined by 
theory and experiment and then are compared with experimental results at low and high energies. 
It is shown that the minimum in the low region is a consequence of a combined effect from the 
vanishing of the vector form factors, the muon mass and Pauli blocking. Several improvements for 
the future investigations are suggested. 
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I. INTRODUCTION 

One pion production is a process, which, along with 
the quasi-elastic scattering and DIS, contribute to the 
total cross section of neutrino interactions with nuclei. 
If the pion is absorbed in the nuclear medium of the 
target, this process constitutes the major background to 
the quasi-elastic process. It is known that the one pion 
production proceeds through resonance production, the 
leading contribution coming from A"*""*" resonance. In 
the coming era of long baseline neutrino experiments the 
cross section of this process needs to be calculated with 
high accuracy. At the same time one must identify and 
separate the coherent component. 

Theoretical and experimental study of one-pion neu- 
trino production was performed in 70's. A comprehen- 
sive experimental study, including the Q^— dependence of 
the differential neutrino cross section, was made in two 
experiments at ANL and BNL. Below we represent the 
data from these two and other experiments and fit them 
with theoretical formulas. 

Theoretically the model of Rein and Segal is of- 
ten used to estimate the neutrino production of the res- 
onances. This model is based on the quark harmonic 
oscillator model for the form factors developed by Feyn- 
man, Kislinger and Randal Another approach is to 
parametrize the neutrino-nucleon-resonance vertex with 
phenomenological form factors. The cross section is usu- 
ally expressed in terms of helicityamplitudes 0, Q with 
detailed formulas given in Ref. [5|. After the discovery 
of neutrino oscillations the production of resonances by 
muon- amd tau- neutrinos was studied again (aSSSIlfl 
as a means for extracting oscillation parameters. The 
calculations have been done, neglecting the muon mass, 
which is a valid approximation for >> to^ . 

Since the time of the ANL and BNL experiments, it 
is known that there is a difference between the data 
and theoretical predictions in the region of small 
(Q2 < 0.1 GeV^). Nowadays it appears that the same 



problem revealed itself in new experiments (K2K and 
MiniBooNE). In this region of the results can be in- 
fluenced by the muon mass, Pauli blocking and coherent 
pion production. In this article we take into account 
the nonzero muon mass and Pauli blocking and calculate 
the cross sections independently. We decided to calcu- 
late the cross section, making use of the phenomenolog- 
ical form factors. Numerical and analytical comparisons 
show, that our results agree with the standard formulas. 

The formulas we present here follow closely the nota- 
tion from deep inelastic scattering where the cross section 
is given in terms of structure fuctions Wi((3^, v), with the 
leptonic variables occuring in multiplicative factors. The 
mass of the muon occurs in the multiplicative factors and 
also enters indirectly when we define the accessible region 
of phase space. We compare our results with the produc- 
tion of A'^^ resonance, where only the amplitude with 
isospin 3/2 contributes. 

The plan of the paper is as follows. In section UTI we 
collect general formulas for the production of resonances 
with J = 3/2. These formulas, together with those in 
Appendix, enable the reader to write a program and pro- 
duce the cross section. We give values for the coupling 
constants and form factors wich are frequently used. In 
section Hill the production of the P33(1232) resonace is 
compared to available data, including muon mass effects. 
Finally, in sections II VI and Ivl wc highlight special proper- 
ties and point out interesting features to be investigated 
in the experiments. 



II. RESONANCE PRODUCTION 

In this article we discuss experiments in which the re- 
action 



(III) 



is studied. We adopt standard kinematics with the defl- 
nitions 



q = k-k', Q^ = -q^, 



= p'^ 



* E-mail 
t E-mail 



olaIakul@zylon.physik.imi-dortmimd.de| 



paschos@physik.uni-dortmund.der 



and compute the cross section ^Qf^^/ ■ The mass of the 
resonance is not restricted to a specific value but allowed 
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to vary within an interval proportional to the width. 
Consequently we let W to vary and write the cross sec- 
tion with formulas analogous to deep inelastic scattering. 
The cross section is now written as 



da 



dVLdE' 



^ cos'' Oc^L^^W^" (II.2) 



with TTiN the mass of the nucleon in the target, M/j the 
mass of the resonance and the leptonic tensor 



and gives the one-dimentional 5— function. Sometimes it 
is convenient to use other variables for resonance produc- 
tion 
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Since the A resonance has an observable width, the 
(5— function should be replaced by its resonance repre- 
sentation 



Li_c,, = rr[7p(l - 75) /zji, /c'] (II. 3) 



The hadronic tensor is defined as 
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where the sum implies a sum over the A polarization 
states and an averaging over the spins of the target. The 
integration over phase space of the A was carried out 
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It is known that resonance production dominates neu- 
trino reactions in the few GeV energy region. The for- 
malism we present in this section is general and holds 
for various resonances. Later on, when we relate the 
structure functions to the form factors, we specialize to 
distinct final states. 

The hadronic matrix element differs from resonance to 
resonance and contains vector and axial form factors. A 
convenient parametrization for the A^'^ resonance is the 
following 



{/\++\r\p) ^^/ii;x{p')d^''u{p) with 
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In the square of the matrix element also appears the Rarita-Schwinger projection operator 
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With these preliminaries the hadronic tensor takes the 
form 



2 2777 JV 



^ ^ - " ^'"'^ -Ml) 

(11.9) 

with (d)'"^ — 7o(c?^)^'^7o and then parametrized accord- 
ing to (|II.4|I . This way we define the relative normaliza- 
tion between the structure functions and the form factors. 
The factor 3 comes from the isospin coefficient for A and 



the 1/2 from the averaging over the initial spins of the 
target. The results of the calculation are summarized in 
Appendix A. 

The remaining problem consists in writing the cross 
section in terms of form factors and specifying their nu- 
merical strength and Q^— dependence. The cross section 
assumes the standard form which includes now the mass 
of the muon. 
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The dependence on the muon mass agrees with the 
one in Ref. j^J. The structure functions are also ex- 
pressed in terms of the form factors. This is straightfor- 
ward and for P33(1232) resonance leads to the relations 
given in Appendix. For final states with opposite parity, 
like 1)13(1520), the 75 matrices in the current-nucleon- 
resonance vertex will appear as multiplicative factors to 
the axial (and not the vector) form factors. The effect 
of this change to Eas. (|A.l() - (|A.6|l is the replacement of 
ruNMR by -mNMn. 

The determination of the form factors follows from gen- 
eral principles and experimental results. We begin with 
the vector form factors. The conserved vector current 
hypothesis gives the relation Cg' = 0. The remaining 
form factors also occur in electroproduction where it has 
been established that the Mi+ multipole dominates. Re- 
cent data determine the contribution from the electric 
multipole E2 to be ~ —2.5% and from the scalar multi- 
pole ^ —5% T^. We shall assume the dominance of the 
magnetic dipole which gives 



remains to introduce the functional form for the width of 
the resonance 



To 



(21+1) 



(11.13) 



with I = 1 for P33(1232) resonance. A partial width of 
the form 



(11.14) 



was also used 0, , but now Eq. HII.13|I is preferable since 
it is required by the partial wave analysis. 

With these results and those in Appendix A, one has a 
complete set of formulas with which to proceed to anal- 
yse the electromagnetic and weak production of the P33 
resonance. 



III. EXPERIMENTAL RESULTS 
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The numerical value is obtained from the data in electro- 
production after an isospin rotation. Electroproduction 
data lead to a Q^— dependence faster than the dipole 
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with Mv = 0.84 GeV. This functional dependence in- 
dicates that the size of the resonance is larger because 
of the mesonic cloud surrounding the resonance and its 
Fourier transform gives a steeper function of Q^. 

Among the axial form factors the most important con- 
tribution comes from whose numerical value is related 
to the pseudoscalar form factor by PCAC. We shall 
use the values 
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with (^A 



15.3, /tt = 0.97m^ and 
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The value for was found to give a small contribution 



to the cross section and the dipole form factor for 

is again modified. The form factor is set to zero as 

suggested by early 0,0,0 dispersion calculations. It 



A detailed experimental study, including the 
Q^— dependence of the differential neutrino cross 
section, was made in two experiments: using the 
Argonne National Laboratory 12-ft bubble chamber 
fANDjlSj and the Brookhaven National Laboratory 7- ft 
bubble chamber (BNL) [TtL H^ . In both experiments 
the neutrino spectrum was peaked at approximately 
1 GeV. We calculate the cross section for the reaction 
vp fi~ — !■ /i^ p 7r+ weighted over the neutrino 
spectrum. 

The BNL experiment observed a peak in the differen- 
tial cross section da/dQ^ at about = 0.175 GeV^ as 
shown in Fig.^ where experimental data are presented in 
arbitrary units. We computed the Q^— distribution with 
the form factors described above and AIa — 1.05 GeV. 
The results are shown in Fig. ^ The theoretical curve 
da/dQ^ has a peak at = 0.085 GeV^ 

It is evident that the agreement between theory and ex- 
periment at > 0.2 GeV^ is satisfactory for both func- 
tional forms of the resonance width (pT3)l and (pT4)l . 
For the overall scale, we normalize the area under the 
theoretical curve for > 0.2 GeV^ to the correspond- 
ing curve under the data. It is also evident that the muon 
mass brings an additional decrease in the region of small 
where the Pauli suppression is also significant, but the 
data are still slightly lower than the theoretical curve. 

In the ANL experiment the data are with large bins 
of and the maximum of da/dQ^ is at a larger value 
of . The formalism described so far determines the 
cross section including the absolute normalization. For 
Ma 



1.05 GeV and the modified dipole in Ea. (|II.12|) 
we obtain the curve in Fig. 13^, which is above the data. 
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FIG. 1: The cross section da/dQ , calculated for the BNL 
neutrino energy spectrum and compared with the experiment 
for the running width (fTTTsll (Fig. (a)) and lUmt (Fig. (b)). 
The full lines are for the case = 0.105 GeV, the dashed 
lines are for the approximation = 0. 



The integrated cross section in this case at high energies 
approaches 0.7 • 10~^® cm^, which is consistent with the 
experimental data. The discrepancy in "dependence 
can not be resolved by the overall normalization of the 
curve and requires a decrease of Ma = 0.84 GeV in order 
to obtain the curve in Fig. The two curves are with- 
out (dotted) and with (solid curve) the muon mass. The 
integrated cross section is also decreased approaching at 
high energies a constant value of 0.55 • 10~^* cm^, which 
is also consistent with the data. 

An earlier theoretical analysis ^3 accounts for the 
ANL data by using similar couplings and muon mass 
effects. They include nuclear corrections by using deu- 
terium wave functions and compare the differential cross 
section to the ANL data. Another approach de- 
scribes electron and neutrino scattering on various nuclei 
in terms of a scaling law abstracted from data and the 
authors present several distributions. A direct compari- 



son with our results is not available and perhaps difficult 
because of the different methods. 

Another way to reach an agreement with the data is 
to replace the dependence (Iii.l2|) with a steeper depen- 
dence, for example 
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+ 

with Ma = 1.05 GeV or 

\0) 



2^2 l + 2Q2/^^r 



Ct{Q^) 



with Ma = 0.95 GeV 

The theoretical formalism is very close to explaining 
the experimental data. There are differences in the BNL 
and ANL data which we can not understand and which 
must be resolved by future experiments. In addition 
there is a sharp decrease at small < 0.2 GeV^ which 
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FIG. 2: The cross section da/dQ , calculated for the ANL 
neutrino energy distribution. The full lines are for = 
0.105 GeV, the dashed lines are for the approximation = 
0. 
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appears to be persistent. One part of the decrease comes 
from the Pauh suppression, which is smaU for deuterium 
and another part from terms depending on the muon 
mass which are important for the low energy of the neu- 
trino beam. 

One can easily see from Fig. |2 that for low energies, 
i.e. neutrino energy i? ~ 1 GeV, taking into account the 
nonzero muon mass reduces the cross section at small 
by approximately 20%. The physical origin of this 
reduction is as follows. Firstly, the double differential 
cross section da/dQ^dW changes mainly due to the con- 
tribution from Ws and structure functions. Secondly, 
for each we must integrate this cross section over W. 
The lower limit of integration W- {Q^) — rriN + is in- 
dependent of and the muon mass. The upper limit of 
integration, however, depends on and the muon mass 
and is given by 



WW) 




1 



+ ^rny 



where s — 2mjv£' + ti^, a± — 1± m'j^/s. 

The integration limits are shown in Fig. O for £' = 
1 GeV and E — 7 GeV. One could easily notice that 
taking into account the muon mass noticeably decreases 
W+ and implies the reduction of da/dQ^. This effect 
diminishes as the neutrino energy increases. In the rest 
of the section the muon mass is taken into account, but 
not discussed any more. 

For the sake of completeness we mention experiments 
at high energies. In these cases effects from the muon 
mass are diminished so that the levelling of the da/dQ^ 
distribution should disappear. The SKAT 2^ experi- 
ment had an average neutrino energy oi E = 7 GeV. 
The results are shown in Fig. Altogether with theoretical 
curves with Ma = 1.05 GeV (which we call case (1)) and 
with Ma = 0.84 GeV (case (2)). There are few exper- 
imental points and the error bars are too large to draw 
conclusions. 

In the FNAL 15-ft bubble chamber |23| experiment 
data are available for neutrino energies between 15 and 
40 GeV. At such energies the integrated cross section 
remains constant with high accuracy, so the exact value 
of neutrino energy is not important. The data and the 
theoretical curves are presented in Fig. |31 

From the BEBC experiment, obtained with the CERN 
wide-band beam, two data sets are available: BEBC-86 
with {E) = 24.8 GeV and BEBC-90 (the neu- 
trino flux is given in ^) with {E) = 54 GeV. They are 
shown in Fig. A common property is the disappear- 
ance of the flattening of the cross section at small Q^, as 
expected. 

We conclude, that different experiments, performed 
with the help of bubble chambers in the 80's, show at 
low a slightly lower cross section than theoretically 
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FIG. 3: Upper integration limit W+(Q'^) for E = 1 GeV and 
E = 7 GeV The full lines are for m,^ = 0.105 GeV, the dashed 
lines are for the approximation = 0. 



predicted. The experiments described so far are not de- 
tailed enough to allow separation of the form factors and 
a unique determination of their dependence. 

Two new experiments K2K and MiniBooNE will be 
delivering results. They are both at low neutrino energies 
where the muon mass effects should be important. With 
the neutrino spectra from and (2^ we predict the 
Q^— distributions shown in Figs. [7| and |H1 using the axial 

These experiments use medium 



form factor in Eq. (|II.12|) 
or heavy nuclei as targets and nuclear corrections must 
be applied, which were left out in our curves (only Pauli 
blocking is included). 



IV. SPECIAL PROPERTIES 

It is evident from our presentation that the cross sec- 
tion in the A resonance region has several important fea- 
tures still to be investigated. One of them deals with 
the structure of the form factors, especially the axial 
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FIG. 4: The cross section da/dQ2, calculated for the SKAT 
experiment with (E) = 7 GeV for the behavior of the form 
factors in cases (1) and (2). 
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form factors. We showed that the dominant contribu- 
tion comes from and Cq. Closely related is the 
(5^— dependence of the process, especially in the small 
region, where various effects are present: 

(i) the mass of the muon is important at low energies, 

(ii) nuclear corrections are also important, like the 
Pauli factor. Up to now there is no sign of a possible 
distortion of the angular distribution due to charge ex- 
change effects, and 

(iii) effects from the mass of the muon should become 
more evident in the MiniBooNE and K2K experiments. 

In order to see the relative importance of the form 
factors we computed in Fig. |^ the various contribu- 
tions for = 1 GcV. The dominance of and 
is evident with the other terms contributing less than 
2 • 10-4° cmV GeV^ As ^ only the contribu- 
tion to da/dQ^ is dominant and remains large for very 
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FIG. 6: The cross section da/dQ2 in experiments BEBC-86 
and BEBC-90 for the behavior of the form factors in cases (1) 
and (2). 



to zero. This is caused by the vanishing of the phase 
space, as we mentioned earlier. Fig.l^D shows the contri- 
bution of the smaller form factors and the sum of them 
is negative with a negative valley at ~ 0.15 GeV^. 

A similar study for the structure functions is shown in 
Fig. El for two values of the neutrino energy (i?^ = 1.0 
and 2.006^). We note that the terms from and 
are negative and W5 contributes to the sharp decrease of 
the cross section at small values of Q^. 

Comparison between experiments with neutrino and 
antineutrino beams will be interesting. The cross sections 
for — !■ /i-pTT"*" and vn iX^mT~ are related to each 
other by changing the sign of the Wa term, i.e. the axial- 
vector interference term. It follows from Fig. that the 
Q^— distribution will be, at low energies, very different. 



small values of Q ~ 0.01 GeV and then turns sharply 



Similarly the difference a^p — crp„ is large at low energies 
and becomes smaller at high energies. This property fol- 
lows from the fact, that the structure of the form factors 
limits the dominant contribution to the low region of . 
Furthermore, the value oi q ■ p = vrriN remains small 
in the resonance region. These two properties together 
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FIG. 7: The cross section da/dQ2, predicted for MiniBOONe 
experiment, W < 1.4 (solid line) and W < 1.6 (dashed line). 
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FIG. 8: The cross section da/dQ2, predicted for K2K exper- 
iment, W < 1.4 (solid line) and W < 1.6 (dashed line). 
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SUMMARY AND CONCLUSIONS 



imply that the contribution from increases quadrat- 
ically with E^^ while the vector-axial interference terms 
grow linearly. To illustrate this property, we define ratio 



da,. 



dan 



R = 



dQ2 dQ^ 



da,. 



dan, 



(IV.l) 



dQ2 rfgs 



and plot it in Fig. ^2 as a function of for three en- 
ergies. The theoretical curves terminate at values of 
when the phase space is kinematically not allowed. 

The channel 9n fi^n'K~ was observed in the BEBC- 
90 experiment for an average antineutrino energy — 
40 GeV. Combining neutrino and antineutrino reactions 
we plot in Fig. ^] the ration i? as a function of Q^. For 
the high energies under consideration and the experimen- 
tally accessible region of < 3.5 GeV this ratio grows 
slowly from 0.0 to 0.15. 



In this article we have calculated general formulas for 
the production of J — 3/2 resonances by neutrinos, in- 
cluding the muon mass. The formalism is applicable to 
the ^33(1232), Z?i3(1520) resonances provided the form 
factors are available. In this article we analyse the pro- 
duction of P33 resonance and compare it with data that 
are available. The analysis of Di3(1520) and the other 
J = 1/2 resonances is left for a future publication. 

Combining results from electroproduction and previ- 
ous analysis of neutrino production we find again that 
the dominant contribution comes from the (Q^) and 
C^{Q^) form factors. Their dependence in is faster 
than the dipole form factors. 

A peculiar feature of the low energy data is a decrease 
of the differential cross section, da/dQ^, as 0. We 

presented an analysis including the mass of the muon 
and Pauli blocking, both of which bring better agreement 
with the data. We have demonstrated that the mass 
of the muon restricts the phase space for the process in 
the low region. The effects should be observable in 
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FIG. 10: Contribution of structure functions to the differen- 
tial cross section. 



the MiniBooNE and K2K experiments. In the small 
region coherent scattering may also be present and should 
be established as a sharp peak in the da/dQ'^dt versus t 
distribution, with the four-momentum-transfer squared 
given by 




-(Y.{e..-A)] (V.l) 



This formula is based on zero energy transfer to the nu- 
cleus but includes the muon mass |27l . Then the in- 
coherent sum of the two effects must reproduce the data. 
It is still interesting to analyse other final state chan- 



nels pir" 



. . as well as reactions with antineutrino 



beams. In the low energy experiments, W3(Q^), which 
distinguishes neutrinos from antineutrinos, is large as it is 
shown in Fig. 1101 For completeness we also presented re- 
sults at higher energies where general trends are already 



apparent. For instance, the mass of the muon becomes 
less important at higher energies and the contribution 
from WsiQ^) is smaller. 

A similar analysis is possible for the 1)13(1520) and 
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FIG. 11: Predictions for R for E ^ 1, 3, 5 GeV. 
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FIG. 12: Experiment BEBC-90: predictions for R. 



J = 1/2 resonances Pii(1440), 5'ii(1535), . . . , where in- 
formation for the vector form factors is now available 
[2^ Isnl l from Jefferson Lab. experiments. We plan to 
include this in the second article of this series. 
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APPENDIX A: STRUCTURE OF THE HADRONIC TENSOR 

As we have mentioned, the hadronic tensor is parametrised in the form ljll.4p . The functions Wi, . . . , Wg have been 
calculated from Ea. l|II.9p and led to 



TUN 



with the Vi {Q^ , f) being the following 



-T = [iq-p-QYiQ-P + rn%) + MU{q-p)^ + Q^m% + Q^mNMR)] 



m% 3 



'at 

cYcY 2 



2 

+ 3 



3 iq-p~Q^)[iq-p-Q^)iq-p + m%~2mNMR) + Mlq-p\ 

[q-p + + ttinMr] 



(A.l) 



y - (^-3 ) 



(CD' 2, 



'AT 



(A.2) 



3 



(~iV f^A 

' '''-{q-p-Q')-CYC^m^ 



TUN 



+ 3(9 -p-Q') 



[2M2 + 2mNMR + Q'^-q-p\ 



/~iv r^A 

'^{q-p-Q')^CXGi 



(A.3) 



These are the important form factors for most of the kinematic region. As mentioned already, there are two 
additional form factors, whose contribution to the cross section is proportional to the square of the muon mass. 



2(G 



V\2 



V - ^^^3)^ [(2? ■p~Q^)iq-p + m%)-Mlim% + tun Mr)] + -'-A^{2q ■ p - Q^) [q ■ p + m% - tunMr] 

o oiVl ^ o Tfl 



V nv 



+ 



2 C^C, 
3Mr ruN 



[{2q-p-Q^){q-p + m%- 2mNMR) + q ■ pM^ 



iCir^+^^i2q-p-Q 



M'r 



+2CtCi - 2^^q ■ p - 2^-^{Ml + Q'-q-p) 



'■AT 



Ml 



[q-p + m\ + vimMr] 



(A.4) 
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2 (^3^)^ r 



3Mfl m,N 



Mi 



2 {CX? r ^ 



'AT 



(A.5) 



(-AT 



O4 Og 



m 



N 



Ml 



{q-p- 



\q-p 



niNMii] 



Va = 



(A.6) 



Notice, as it is expected, for the contribution of the 
vector form factors the equahties W5 = W2 • (<Z ■ p)/Q^ 
and W4 = W2 • (g • p)^/Q^ - Wim%/Q'^ are satisfied. 

In terms of the invariant variables, and W , the 
scalar products of the 4- vectors are: 



k-p — niNE, q-p = m,Ni', v 



2m 



N 



the Ref.[i. This can be easily done for a specific case 
0, — > 0, when only contribute to the cross 
section. After the integration over W with the help of 
the delta-function we obtain 



da 



ZTT 



1 - 



An 



m 



N 



2mNE 



2Ml 



p-p ^-{W^ + + m^^), q-p'^mNU-Q\ 
Eq. pi.5|l must be compared with the known one from 



that identically coincides with the result, obtained in a 
similar way from 
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